Correlated random band matrices: localization-delocalization transitions 
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We study the statistics of eigenvectors in correlated random band matrix models. These models 
are characterized by two parameters, the band width B(N) of a Hermitian N x N matrix and the 
correlation parameter C(N) describing correlations of matrix elements along diagonal lines. The 
correlated band matrices show a much richer phenomenology than models without correlation as soon 
as the correlation parameter scales sufficiently fast with matrix size. In particular, for B(N) ~ \/~N 
and C(N) ~ \/jV, the model shows a localization-delocalization transition of the quantum Hall type. 



PACS numbers: 02.50.-r Probability theory, stochastic processes, and statistics; 73.23.-b Mesoscopic 
systems; 71.30.+h Metal-insulator transitions and other electronic transitions 



mm m • u 

pi i i m i 

□ ■BW8M 

■■ JEMS 
■■□■■□■□■□BSD 

n , ■ : _ m m ■■ 

JDIIIGIIIIUII 

■ □■□□BBBBtBHDnra 
EE Till ■ 

mum: : iiim ^ 
■ mam jhhbh 
■■ill " r f 
'■ 



□ □□□ 

□ CHD 

□ □□ 

□ □□□ 

□ □□□ 

□ □□□ 

□ □□□ 

□ □□□ 

□ □□□ 

□ □□□ 

□ □□□ 

□ □□□ 



□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□□ 
□□□□□□□□ 

□□□□□□HDD 

□□□□□□□»□ 
□□□□□□□(JO 



■ S[ ][ ] 
' J ■ 

□hlm 

■ 

□□□□ 
□□□□ 
□□□ 

□□ 

,. Ml 
□ □□ 



[!□□□ □□! ;' 

: :: : 
: :: :: :: : 

□□□□□■Hiiri 

uuuuuumsn 
□□□□□□□an 



N=13 
of 

3x3x3 



(a) 



(b) 



FIG. 1. Visualization of typical random matrices obeying 
Wigner-Dyson statistics (a), Poisson statistics (b). The in- 
tensity represents the absolute value of matrix elements. 



FIG. 2. Visualization of a typical random matrix obeying 
the statistics of a 3D Anderson model. The intensity repre- 
sents the absolute value of matrix elements. 



I. INTRODUCTION 

In the theory of random Hermitian matrices M two ro- 
bust types of statistics are found in the limit of infi- 
nite matrix size (denoted here as 'thermodynamic limit'). 
First, the Wigner-Dyson statistics describing systems 
that become ergodic in the thermodynamic limit and 
have an incompressible, correlated spectrum and Gaus- 
sian distributed, uncorrelated amplitudes of the corre- 
sponding eigenstates (see Fig. la). Since we do not con- 
sider further symmetry constraints we focus on the ma- 
trix ensembles denoted as class A in the classification 
of [||. A convenient representative of its ergodic lim- 
iting ensemble is given by the Gaussian unitary matrix 
ensemble GUE. The second robust statistics is the Pois- 
son statistics with eigenstates, localized on certain ba- 
sis states (sites), and with an compressible, uncorrelated 
spectrum (see Fig. lb). Real complex quantum systems, 



represented by random Hermitian matrices, can show 
a crossover between Wigner-Dyson and Poisson statis- 
tics or, in some cases, a true quantum phase transition 
with novel 'critical' statistics. A well known example 
is the 3D Anderson model (see Fig. 2) describing the 
motion of independent electrons on a 3D lattice with 
random uncorrelated on-site disorder. Below a certain 
critical value of disorder, in the thermodynamic limit, 
all states at the energy band center are infinitely ex- 
tended (delocalizcd) in space, while for larger disorder 
all states are spatially localized. Instead of changing the 
disorder, one can change the energy within the energy 
band, keeping the disorder fixed at low values. Again, 
a transition from localization (band tails) to derealiza- 
tion (band center) occurs. It is worth mentioning, that 
the average density of states (DOS) is non-critical, i.e. it 
stays smooth across the localization-delocalization (LD) 
transition. Although these statements are substantiated 



1 



by analytical as well as numerical work (for reviews see 
HQ), the special structure of this matrix ensemble (com- 
posed of a sparse, but deterministic matrix and a ran- 
dom diagonal matrix) has prohibited, so far, a rigorous 
proof of these statements. Another well known system 
with a transition from localized to critical states is the 
two-dimensional (2D) quantum Hall system (for reviews 
see ]5],[|) which we will describe briefly later. Further- 
more, several matrix ensembles modeling the motion of 
2D disordered electrons undergoing (time-reversal sym- 
metric) spin-orbit interactions are known to display a LD 
transition (see e.g. J7J). In all of these realistic matrix en- 
sembles the statistics at criticality represents an unstable 
fixed point under increasing system size (i.e. matrix di- 
mension), which means that any slight shift away from 
the critical value of the energy, say, will drive the system 
into one of the stable matrix ensembles, Wigner-Dyson 
for the delocalized states and Poisson for localized states. 
The critical ensembles are characterized by correlated 
spectra, but with a finite compressibility. Furthermore, 
critical eigenstates are multifractal and the multifractal 
exponents are related to the compressibility of the spec- 
trum (for a review see Q). 

It is desirable to study matrix ensembles with simple con- 
struction rules and to ask for necessary ingredients in or- 
der to have a LD transition. Also in quantum chaos the 
interest in crossover ensembles has grown B. In that 
context the Rosenzweig-Porter model Q was studied as 
a toy-model for the crossover. It is defined as a simple su- 
perposition of a Poissonian and a Wigner-Dyson matrix. 
It has been shown rigorously that, by choosing the super- 
position in an appropriate way, novel critical ensembles 
emerge, but the spectral compressibility is identical to 
the Poisson ensemble and states are not multifractal (see 
p0| and references therein). Another well studied matrix 
ensemble is that of random band matrices (RBM) with 
uncorrelated elements. The band width B describes the 
number of diagonals with non-vanishing elements. For 
B ~ N s , with s > 1/2, one recovers the Wigner-Dyson 
statistics. Such band matrix models have been discussed 
in the context of the 'quantum kicked rotor' problem Jll]] 
and have been studied extensively in a series o f pape rs by 
Mirlin, Fyodorov and others (for review see [[i2||l5f |). It 
turned out that, in particular for B 3> 1, all states are lo- 
calized with a localization length (in index space) £ ~ B 2 . 
For fixed B one has therefore a crossover from Wigner- 
Dyson to Poisson statistics as N is taken from values 
much smaller than B 2 to values much larger than B 2 , 
and B 2 /N is the relevant parameter for a scaling analy- 
sis of data. Superpositions of such random band matri- 
ces with random diagonal matrices have been studied in 
the context of the 'two-interacting particle' problem (see 
e.g. |l4|,|l5|] ) , however these ensembles do not show novel 
critical behavior as compared to the Rosenzweig-Porter 
model. 

In fact, only few simply designed matrix ensembles are 
known to become critical with multifractal critical states 
(see Q,^), for example 'power law' band matrix en- 



sembles, where the strength of (uncorrelated) matrix el- 
ements falls off in a power law fashion in the direction 
perpendicular to the central diagonal. The critical cases 
occur for the power law behavior ~ a; -1 of the typical ab- 
solute values of matrix elements |16|Jls[] . It is, however, 
important to notice a significant difference to realistic 
critical ensembles: there is no LD transition within the 
spectrum; iff parameters are fixed to critical values all 
states are critical. 

In this paper we study correlated random band matrix 
(CRBM) ensembles and, with the assistance of numeri- 
cal calculations, argue that these ensembles can lead to 
a LD transition within the spectrum. A new parameter 
C(N) ~ iV* describing the correlation of certain matrix 
elements is introduced for random band matrices. For 
B(N) ~ \/~N states are localized outside of the energy 
band center and a LD transition at the energy band cen- 
ter occurs provided 1/2 < t < 1. 

A major motivation for studying these ensemble origi- 
nates from the theory of the integer quantum Hall effect 
The plateau to plateau transition in the quantum 
Hall effect can be captured in models of non-interacting 
2D electrons in a strong magnetic field and a random 
potential, referred to as quantum Hall system (QHS). In 
the one-band Landau representation the Hamiltonian is 
represented as a random matrix with two characteristic 
features, (i) The matrix elements decay perpendicular 
to the main diagonal in a Gaussian way. (ii) No correla- 
tions exist between elements on distinct 'nebendiagonal' 
lines, but Gaussian correlations exist along each of the 
nebendiagonals. These features led to the introduction 
of the 'random Landau model' (RLM) to study critical 
properties of QHSs (see e.g. ||). The original purpose 
of constructing the RLM was to avoid explicite calcula- 
tions of matrix elements starting from a randomly chosen 
disorder potential and to directly generate the matrix ele- 
ments as random numbers that fulfill the statistical prop- 
erties (i) and (ii). As will be explained in more details be- 
low, the corresponding CRBM model simplifies the RLM 
further, in as much as a sharp band width is introduced 
and correlations along nebendiagonals are idealized and 
cutoff after a finite length. 

Recently, matrix ensembles with correlated matrix ele- 
ments attracted some interest [|l9] in the context of the 
metal-insulator experiments in 2D ]2C|], for which strong 
Coulomb interaction is believed to be a necessary ingre- 
dient. It is very interesting that in |2|] ID models with 
correlated disorder potentials could, to a large extent, be 
solved analytically and that certain correlated disorder 
potentials were shown to cause LD transitions in ID. It 
is not obvious how to extend the method of |^l| to the 
case of CRBM models. Usually, correlations in matrix 
ensembles lead to serious complications in analytical at- 
tacks. For example, in the field theoretic treatment (see 



e.g. 1 22 1) of random matrix ensembles the absence of long 
ranged correlations is essential to find appropriate field 
degrees of freedom that depend smoothly on a single site 
variable. In our CRBM models correlations are intro- 
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duced by constraints (a number of matrix elements are 
taken to be identical). This may help to reduce com- 
plications in constructing a field theoretic approach for 
CRBM models. 

In Sec. 2 we give a detailed definition of the CRBM and 
discuss three alternative interpretations. The investiga- 
tion of the LD transition is carried out in Sec. 3 by a 
multifractal analysis of states for an ensemble that is ex- 
pected to fall into the quantum Hall universality class. 
Our results are in favor of this expectation. The analysis 
is carried over to ensembles where correlations are taken 
to extreme limits in Sec. 4. In Section 5 we present our 
conclusions. 
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II. CORRELATED RANDOM BAND MATRIX 
MODEL 

Let the elements of a N x N Hermitian matrix H be 
written as 

H k i =x M +iy kl for I > k , (1) 
H kk = V2x kk , (2) 

where all non- vanishing real numbers x k i , y k i are taken 
from the same distribution V with vanishing mean and 
finite variance a 2 . We take the symmetric and uniform 
distribution on [-1,1] (a 2 = 1/3). With B, C being two 
integer numbers, called the 'band width' and the 'correla- 
tion parameter', respectively, the correlated band matrix 
ensembles are defined by the following algorithm (I) 
(III) and are visualized in Fig. 3. 

(I) Begin with the main diagonal of H and draw a random 
integer number N\ < C, and the random number x\\ 
from the distribution V . Take the first N\ values on the 
main diagonal equal to \[2x\\. Then, draw xn 1 +i,n 1 +i 
from the distribution V and take successively C elements 
on the main diagonal equal to \/2(xiVi+i,jVi+i)- Now, 
draw arjVi+c+i,JVi+c+i from the distribution V and take 
successively C elements on the main diagonal equal to 
\/2a;jVi+c+i,iVi+c+i- Continue with this procedure until 
the main diagonal is filled up. 

(II) Now consider the next 'nebendiagonal'. Its first ran- 
dom clement H12 = x i2 + iy\2 and a random integer 
number N 2 < C are drawn. Put the first N 2 elements of 
this nebendiagonal equal to Hyi- The next C elements on 
the same nebendiagonal are taken equal to the random 
value -ffw 2 +i,A r 2+2, and so on. 

(III) The procedure terminates after the -Bth nebendiag- 
onal (I — k = B) is filled up. All other matrix elements 
(I — k > B) are set to zero. Finally, Hermiticity is in- 
stalled by taking Hi >k = H% <t . 

For C = 1 the usual band matrix models (band width 
B) of uncorrelated matrix elements are recovered. C 
and B form the relevant parameters of the CRBM, while 
the value of a is not significant it just defines the 
energy units. For finite C > 1 the correlation along 
(neben)diagonals is a triangular function of range C and 



FIG. 3. Visualization of a typical matrix of the correlated 
random band matrix model. The intensity represents the ab- 
solute value of matrix elements. 




FIG. 4. One-dimensional interpretation of correlated band 
matrices. Hopping events for distinct hopping distances are 
uncorrelated. Hopping events for the same hopping distance 
are correlated over typically C/2 nearest neighbors. In this 
example B > 6, C > 3. 

half-width C/2. Thus, C/2 is a typical distance over 
which elements are correlated along (neben)diagonals. 
The spectrum is always distributed in a symmetric way 
around the center E = as a consequence of the symme- 
try of the distribution V . 

In the following, we are going to discuss three possible 
physical interpretations. The most obvious interpreta- 
tion relies on the site representation | I) — (0, 0, . . . , I = 
1,0,..., 0). In this representation, H describes hopping 
of particles on a ID chain of length L — N (lattice spac- 
ing = 1) with a maximum distance of hopping equal to 
B. The average hopping probability in one instant of 
time t is {% = 1, (. . .) denotes the ensemble average) 

(| (k{t)\l{t + 1)) | 2 ) =a 2 for < \l - k\ < B . (3) 

The correlation between matrix elements means that two 
hopping amplitudes are identical, if the hopping distance 
is equal, and provided the hopping starts at sites the 
distance between which is less than typically C/2 (see 
Fig. 4). 

An alternative interpretation results, when the N sites 
are arranged in a quasi-lD (Q1D) geometry with N c = 
B+l parallel 'channels' in a 'wire' of length L 1 = N/(B + 
1). The lattice spacing along the wire, h, and the unit 
of time arc taken to be 1 . Transitions in one unit of time 
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FIG. 5. Quasi-one-dimensional interpretation of correlated 
random band matrices. Sites are labeled such that for a fixed 
cross section of the wire all channels are filled up with sites, 
before one proceeds to the next cross section. Hopping events 
along the wire and couplings within each cross section are 
uncorrelated, if label distances are distinct. For identical label 
distances they are correlated over typically C/2 sites. In this 
example B — 3, C > 2. 

are possible between channels (coupling) and along the 
wire (hopping). Sites are labeled (see Fig. 5) such that 
hopping is possible at most over one lattice spacing along 
the wire. Again, hopping (coupling) is correlated over 
typically ~ C/2 states, provided the difference between 
labels is identical. A less obvious interpretation of CRBM 
models arises when quantum Hall systems are studied 
in a one band Landau representation. A quantum Hall 
system is described by a one-particle Hamiltonian of 2D 
electrons (charge — e, mass m) in the presence of a strong 
perpendicular magnetic field B and a random potential 
V{x,y). In Landau gauge the Hamiltonian reads 

H = i ( p * + iPy + eBx)2 ^ + v{x ' v) ' (4) 

where (j> x ,p y ) is the canonical momentum with respect to 
the Cartesian coordinates (x,y). For periodic boundary 
conditions in y-direction (length L y ) the kinetic energy 
forms a highly degenerate harmonic oscillator problem 
(p v is conserved) that is diagonalized by Landau states 
p3[ | n, /). Here n = 0, 1, 2, 3, . . . labels the Landau en- 
ergies E n — TuOcin + 1/2) (lo c = eB/m cyclotron en- 
ergy), and I = 0, 1, 2, . . . labels center coordinates of the 
degenerate Landau states. The Landau states are sepa- 
rated into plane waves in y-direction with quantized mo- 
mentum qi and into oscillator wave functions centered at 
Xi = —\ 2 qi, where A = y/h/ (eB) is the characteristic 
'magnetic length'. As long as the typical values of the 
random potential are much smaller than the cyclotron 
energy, one can study the full eigenvalue problem of the 
low lying 'Landau bands' approximately by restricting 
the Hilbert space to separate Landau levels n. In par- 
ticular, for the lowest Landau level the Landau states 
read 

<x ' W) = exp (- p i — 2A3-J ' 

(5) 

A convenient recipe to study finite systems of length L in 
a;-direction is to use 'Landau counting' of states, that is 
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FIG. 6. Interpretation of correlated band matrices as quan- 
tum Hall system. Bare states are interpreted as Landau states 
on a 2D plane. The qualitative properties of Landau states 
are shown: plane waves in y direction situated at center co- 
ordinate Xi, extended in rc-direction over approximately one 
magnetic length A corresponding to N\ neighboring states. 

to take only those Landau states into account for which 
the center coordinates Xi fall into the interval [0, L\. The 
total number of states, for an aspect ratio a = L y /L, is 

N= (a/2ir)(L/\) 2 . (6) 

By shifting the lowest Landau energy to zero, the eigen- 
value problem is defined by the matrix H k i = (k\V\l) 
which reads 

H kl = P -^=- dxV(x;AX)e-&*- x) ' , (7) 

V 7 ™ J 

— oo 

Ly/2 

V{x;AX) = L~ x J dyV{x, y yv AX ' x \ (8) 

where X = (X k + X t )/2 and AX = X k - X t . These ma- 
trix elements form a random N x N matrix, the elements 
of which are composed by a Fourier transformation of the 
random potential in y-direction and a Gaussian weighted 
averaging of the potential over a magnetic length in in- 
direction. Crucial for the structure of the matrix is the 
fact that, within the distance of one magnetic length A, 
a number of N\ ~ N(X/L) different Landau states can 
be situated (see Fig. 6). Thus, for a constant aspect 
ratio, this number increases as N\ ~ s/~N in the thermo- 
dynamic limit. This leads to correlations between ma- 
trix elements along nebendiagonals over a range of N\ 
states. This range increases when the disorder poten- 
tial is correlated in real space over distances exceeding 
the magnetic length. The correlation between matrix el- 
ements on distinct (neben)diagonals is negligible, since 
the correlator contains the superposition of ~ N\ ran- 
dom phase factors. The matrix elements decay perpen- 
dicular to the main diagonal due to the Gaussian decay 
of Landau states. For disorder potentials with a spatial 
correlation length d ^S> A the decay sets in earlier, be- 
cause the Landau states are orthogonal. A quantitative 
analysis shows (see ||) that the Gaussian band- width is 
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and 

the Gaussian correlation length along (neben)diagonals 
is 



C = (3 



(10) 



Here (3 > 1 is a parameter that is controlled by the po- 
tential correlation length d, 



Vd 2 + A 2 
A 



(11) 



For finite and fixed values d, a, in the thermodynamic 
limit, B(N) and C(N) increase as y/N to infinity. A 
matrix for QHSs with the statistical properties described 
above is denoted as 'random Landau matrix' (RLM) (see 

Our CRBM simulates the RLM in as much, as it has the 
same qualitative features of a finite band width B and a 
correlation length, being typically C/2. The quantitative 
differences are that, in the CRBM, the matrices have a 
sharp band width B, instead of a Gaussian band width 

B, and that the correlation function is a triangle with a 
half width C/2, instead of a Gaussian with a half width 

C. In the thermodynamic limit, however, we expect that 
these differences should be insignificant for the statistics 
of eigenvalues and eigenvectors, provided the parameters 
B(N) and C(N) scale in the same way with N as the pa- 
rameters B(N) and C(N), respectively. Note, however, 
that in the RLM the ratio C(N)/B(N) = (3 2 is bounded 
from below by 1, while in the CRBM model we are free 
to choose any value for B(N), C(N). 



III. MULTIFRACTALITY AND SPATIAL 
CORRELATIONS 

To study the LD properties of matrix ensembles nu- 
merically, one can follow a number of different strate- 
gies. The most efficient is to analyze only the eigen- 
value statistics. Although the eigenvalues encode most of 
the relevant information about LD properties, the statis- 
tics of wavefunction amplitudes is more direct. Local- 
ization, derealization and even criticality of states can 
be qualitatively distinguished already by inspection of 
plots of the squared amplitudes of wavefunctions (see e.g. 
Fig. 7). Critical states ?/>(r) are characterized by hav- 
ing a multifractal distribution of its squared amplitudes 
prob(r) = |f/;(r)| 2 . This spectrum becomes independent 
of system size and is universal for all of the critical states 
in the thermodynamic limit (for a review see jgjj]) or, 
more generally, follows a universal distribution ( see pq| ) . 
In particular, the geometric mean is a convenient mea- 
sure of a typical probability and scales as 



where the deviation ao — d > of the fractal dimen- 
sion ao from the Euclidean dimension d signals multi- 
fractality and is the most sensitive critical exponent of 
the LD transition. Although a critical state is extended 
all over the system, it fluctuates strongly and has large re- 
gions of low probability which results in the stronger de- 
cay of typical amplitudes as compared to homogeneously 
extended states. A quantity that is closely related to 
ao is the exponent rj of long ranged spatial correlations 
(prob(r)prob(O)) ~ r _?? pq , p7[ . It fulfills scaling rela- 
tions to the fractal dimension of the second moment of 
prob(r) pH , and also to the compressibility of the eigen- 
value spectrum [^9| (see also [|o|). As a crude estimate 
(based on a log-normal approximation |2$| to the distri- 
bution of prob) one has r\ f=a 2(ao — d). 
In this work we focus on wavefunction statistics and the 
determination of ao- One should, however, be careful 
when drawing conclusions from the calculation of frac- 
tal exponents for finite matrix sizes. Such calculations 
should be assisted by inspection of states, and one should 
study the dependence on the matrix size N. For exam- 
ple, states with localization lengths that are small, but 
not very small compared to system size tend to produce 
larger values of ao because parts of the wave function 
have low amplitudes. This can be seen easily in plots of 
the corresponding squared amplitudes. From the linear 
regression procedure that allows to determine ao one can- 
not distinguish such behavior from true multifractality, 
as long as the system sizes cannot be made much larger 
than the localization length. Furthermore, one should 
distinguish between spatial, energy, and ensemble aver- 
ages. In practice, we first perform the spatial average for 
a fixed wavefunction to determine the exponent ao for 
a finite matrix size N by the box-counting method (see 
e.g. pq|). This can be done for many states in the critical 
region (which has finite width for finite N) and we aver- 
age over these states. Finally, an average over different 
realizations (ensemble average) can be performed. It is 
not obvious, that the order of these different procedures 
commutes since the exponent fluctuates from state to 
state. The question about the self-averaging of the frac- 
tal exponents was recently addressed in p5j and it was 
claimed that exponents follow a universal scale indepen- 
dent distribution function in the thermodynamic limit, 
rather than being self-averaging. We do not investigate 
this question here. For our finite systems the exponents 
are fluctuating anyway and we consider averages (over 
typically 100 states) as explained above. 
For a number of different models of QHSs the exponent 
ao has been determined, e.g ao = 2.28 ± 0.02 in jn). 



prob t = exp (ln(prob(r))) - L~ 



(12) 



Other authors (e.g. |32|) find values between 2.27 ± 0.1 
and ao = 2.29 ± 0.02. The largest system sizes studied 
were about L sa 200A leading to matrix dimensions N of 
about 10 4 . 

For the sake of a direct comparison of wavefunction plots 
between our CRBM models and the RLM we took the 2D 
Landau representation of Eq. (^|) with a Landau counting 
for an aspect ratio a — 1. The eigenvalues and eigenstates 
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FIG. 7. Squared amplitudes of a localized state (a), and a 
critical state (b), in the correlated random band matrix model 
for the standard quantum Hall case. States are represented 
in two-dimensional Landau representation. 



were calculated numerically exploiting the band structure 
of the matrices. To be as close as possible to a real QHS 
with an aspect ratio awl and a short ranged random 
potential, we took B(N) = B(N, a = 1,0 = 1) = ^/N/n 
and C(N) = (l/2)C(N,a = 1,0 = 1) = {l/2)^/WJ^ 
(e.g. the largest matrices had parameters N = 6400, B = 
45 and C = 90). We refer to this choice of parameters as 
the "standard quantum Hall case" 

Our findings for the standard quantum Hall case of the 
CRBM model can be summarized as follows. Almost all 
states are localized (see Fig. 7a) . Only those around the 
energy band center E = are multifractaly extended 
(see Fig. 7b). For finite N there is a small energy band 
of extended states with localization length £ > L. The 
energy band width of extended states, A c , shrinks with 
increasing TV with a critical exponent related to the di- 
vergence of the localization length at the band center. To 
determine this critical exponent it would be sufficient to 
determine A c as a function of system size. However, A c 



can be defined only by considering an appropriate scaling 
variable, e.g the participation ratio, that increases above 
some threshold value when the states become extended. 
In finite systems scaling variables are strongly fluctuat- 
ing and one has to consider distribution functions and/or 
appropriate typical values. We did not try to calculate 
this exponent precisely, but only convinced ourself that 
the typical numbers of clearly extended states was com- 
parable to those of realistic quantum Hall systems with 
B(N) — B(N), C(N) — 2C{N). 

The fluctuation of ao determined by the box-counting 
for individual states can bee seen in Fig. 10. The average 
over 100 different extended states of a N — 6400 standard 
quantum Hall case yields ao = 2.26 ± 0.02. This value is 
close to the value ao = 2.28 ± 0.02 obtained by the same 
ave rag ing procedure for an original quantum Hall system 
in plj] . We therefore conclude, that the CRBM shows 
indeed a LD transition reminiscent of quantum Hall sys- 
tems. It seems likely that, in the thermodynamic limit, 
the critical behavior of CRBMs in the standard quantum 
Hall class is actually identical to that of the quantum Hall 
universality class, because the essential TV dependence of 
the relevant parameters, band width B(N) ~ ^/N and 
correlation parameter C(N) ~ V~N, are identical to those 
of realistic QHSs in Landau representation. 
So far the comparison of multifractal exponents was 
based on the wave function statistics without any refer- 
ence to spatial correlations. Therefore, a more ambitious 
comparison between the CRBM and a QHS concerns the 
critical exponents of spatial correlations and their scaling 
relations to the multifractal exponents. In a multifrac- 
tal state the g-correlator (prob 9 (r)prob 9 (0)) ~ r^ x ^ has 
fractal dimension 



x(q) = 2A(q) - A(2q) 



(13) 



where A(g) are the usual fractal exponents of the q- 
moments (prob 9 (r)) ~ L _A («) (for review see H). 
We find for the spatial correlations of the standard quan- 
tum Hall case the scaling exponents shown in Fig. 8. 
They are compared with the data that follow from the 
spectrum A(g) and Eq. (13). The spectrum A(g) was 
calculated by the box counting-method. They satisfy, 
within the numerical uncertainties, the scaling relation 
Eq. (|lj). Furthermore, the exponents are close to their 
values for critical states in realistic quantum Hall sys- 
tems 1 31 , e.g. the exponent of 'anomalous diffusion' |27|j 
isa;(l) =0.4 ±0.1. 

We summarize this section by stating that our multi- 
fractal analysis gives a number of clear indications that 
correlated random band matrices, in the "standard quan- 
tum Hall" case, are true representatives of the quantum 
Hall universality class. 

We close this section with an observation not directly rel- 
evant for the questions addressed in this paper, but which 
may be relevant for those readers that like to perform nu- 
merical calculations with the CRBM model. In realistic 
QHSs with an aspect ratio a ^ 1 one observes that a 
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FIG. 8. Numerically obtained correlation exponents z(q) 
(symbols) in comparison with the data following from the 
scaling relation (jl3|) (line). Error bars are shown for two 
representative values. 

number of the multifractal states tend to have an orien- 
tation along the direction of smaller width. This orienta- 
tion effect has no influence on the asymptotic statistical 
properties of the wave function amplitudes as N scales to 
oo. Only corrections to scaling due to finite system sizes 
can be different for different aspect ratios. However, an 
aspect ratio a ^ 1 will influence scaling variables like the 
Thouless sensitivity gxh = Se/A. It measures the change 
in energy Se due to a change from periodic to antiperiodic 
boundary conditions in a given direction, relative to the 
mean level spacing A. It is exponentially small for local- 
ized states and typically of order 1 for critically extended 
states. We calculated this quantity for a realistic quan- 
tum Hall system. It has strong mesoscopic fluctuations 
(variance ~ mean) and found that the unique maximum 
of its typical values at the band center, for a = 1, splits 
up into two maxima, for a / 1, symmetric around the 
band center. This behavior is related to the fact that 
those wave functions that start to be extended in the 
direction of smaller width are more sensitive to changes 
in the boundary condition than those that have already 
huge localization lengths and are uniformly extended in 
both directions. In the CRBM model we observe a sim- 
ilar phenomenon. For the standard quantum Hall case 
we actually found a tendency for an orientation into the 
y-direction (with Landau states as defined in Eq. (0) and 
Landau counting for a = 1)). This behavior changes to 
an orientation in the opposite direction under increasing 
C(N) by a factor of 0(1), keeping B(N) fixed. In con- 
trast to a realistic QHS where we can calculate B, C for a 
truly symmetric situation, a — 1, in the CRBM model we 
do not know a priori if the choice C/2 = B is appropriate 



to simulate a truly symmetric situation, a = 1. Because 
the identification of the half-width value of a triangular 
correlation function with the half-width of a Gaussian is 
not strict, taking a factor of order unity between them 
is equally well justified. The same ambiguity is present 
in the identification of the band width B with B. Any 
change in B, C by a factor of order unity can therefore 
lead to the orientation effect. As in realistic quantum 
Hall systems we also found the splitting of the maximum 
in Thouless sensitivities in the standard quantum Hall 
case. This splitting effect may be further investigated. 

IV. TUNING THE CORRELATION PARAMETER 

C{N) 

When the correlation parameter is fixed to a constant 
C(N) = Cq the CRBM will behave like an ordinary ran- 
dom band matrix when A' > Cp. This case is very well 
understood (see Jl2|,|l3|]) and one knows that a crossover 
from localization to Wigner-Dyson delocalization takes 
place as the band width B is varied. The localization 
length in ID interpretation is £id = cB 2 (c a constant 
of order unity). Except for the far energy band tails, 
where states are stronger localized, the localization be- 
havior is almost uniform over the energy band. It is worth 
noting that the amplitudes of wavefunctions within the 
central region (where the amplitudes are not exponen- 
tially small) are strongly fluctuating, but they are not 
multifractal in the limit of large B(N). The entire dis- 
tribution of amplitudes is, asymptotically in B, fixed by 
the value of the ratio g = £ 1D /A = cB 2 (N)/N. This 
ratio is the relevant scaling parameter and is denoted as 
" conductance" . As we have seen in the previous section, 
the behavior of CRBM changes drastically when the cor- 
relation parameter increases sufficiently fast with N. In 
the standard quantum Hall LD transition takes 

place within the energy band. To get more insight into 
the role of the correlation parameter we therefore studied 
two extreme cases: (A) C(N) = 1, and (B) C(N) = N. 
In both cases we kept B(N) ~ yN as in the standard 
quantum Hall case. 

Situation (A) corresponds to the usual uncorrelated ran- 
dom band matrix models with a large localization length 
of £id ~ B 2 ~ L in ID-interpretation (£qid ~ B ~ L' 
in Q ID-interpretation) and a constant 'conductance' 
of order 1, g = £m/L — £qid/L' = const. This 
model has no interpretation as a QHS, since the ratio 
C(N)/B(N) ~ N' 1 / 2 < 1. For better comparison we 
used the same Landau representation as before and per- 
formed a multifractal analysis of the extended states by 
the same box-counting method as in the standard quan- 
tum Hall case. 

Our findings in situation (A) can be summarized as fol- 
lows. All states behave similar within the band (except 
for those in the far tails) The states are not uniformly 
extended, but are confined to strips with a width of 



7 



0.0006 
0.0005 
0.0004 
0.0003 
0.0002 
0.0001 





20 40 60 80 100 



FIG. 9. Squared amplitudes of a typical state in the ran- 
dom band matrix model for situation (A) with vanishing cor- 
relation between matrix elements. The state is represented in 
the two-dimensional Landau representation. 



about half the system size (see Fig. 9). Within that strip 
the states are extended and they fluctuate strongly in a 
'grassy' way. They do not show the self-similar regions of 
low amplitude like typical multifractal states. This be- 
havior is compatible with the ID (or quasi-lD) interpre- 
tation of the uncorrelated random band matrix with a lo- 
calization length of the order of L (£') and a conductance 
of order unity. We calculated, for N = 6400, B — 45, the 

\1D I (J— 11 

exponent q;q j w 2.14 (see Fig. 10). This value must 

be taken with care, as the states were not extended over 
the full system. They are localized to an area of about 
half the system size. Thus, the regions of exponentially 
small amplitudes outside the localization center lead to 
values a > 2. Taking amplitudes from only the local- 
ization center reduces the average of ao, but fluctuations 
from state to state are strong. We therefore expect that 
ao, measured in the localization center, will slowly con- 
verge to ao = 2 as B(N) increases further with N. 
Situation (B) deviates strongly from the usual uncorre- 
lated random band matrix models. The elements on each 
(neben)diagonal are constant, but uncorrelated for dis- 
tinct (neben)diagonals. We denote them by h m where 
m=0 labels the main diagonal and positive (negative) m 
label the upper right (lower left) lying nebendiagonals. 
Had we set N = oo in the first place, we could solve the 
eigenvalue problem by Fourier transformation. Each hop- 
ping event over a fixed distance would be translational 
invariant. Thus, the eigenstates, for N = oo, are plane 
waves ip q (l) = e lql where q is a quantum wave number 
that can take any real value. The corresponding eigen- 
value is 



E q = 



E >>« 



= h + 2di 



B 

E 



(14) 



In Landau representation the plane waves ip q (l) trans- 
form into wave functions ip q {x,y) that are plane waves 



FIG. 10. Multifractal exponents for 100 extended states 
of correlated random band matrix models in the standard 
quantum Hall case (middle), in situation (A) (bottom) and 
in situation (B) (top). The straight line corresponds to the 
average in the standard quantum Hall case. 



in x direction, centered at a center coordinate Y q = —\q, 
and have a width of a magnetic length in y direction. 
For any finite N, however, such solution is not possi- 
ble, unless periodic boundary conditions are implemented 
in the site representation. To implement them into our 
band matrix models we have to add ~ B 2 matrix ele- 
ments in the upper right (and lower left) corners of the 
matrix. This would violate the band structure. We see 
that, for any finite N, the correlated band-matrix brakes 
the translational invariance of hopping events, and it is 
not obvious that the states restore this symmetry when 
N goes to infinity. Actually, our finite N results indicate 
that the states will not be plane waves in the center of the 
band (see Fig. 11). Furthermore, a simple perturbative 
treatment shows that the omission of the ~ B 2 elements 
in the corners cannot be neglected in the limit N — ► oo. 
The CRBM in situation (B) does also allow for an inter- 
pretation as a quantum Hall system, since C(N)/B(N) ~ 
V~N > 1. As follows from equations (§- 0) the potential 
correlation length d ~ N 1 / 4 and the aspect ratio is large, 
a ~ \f~N . This translates to the scaling with system size 
L as 



d ~ L ~ 7V 1/4 , L y =aL~L 3 . 



(15) 



The CRBM in situation (B), thus represents a long quan- 
tum Hall strip where L y /L ~ (L/A) 2 and the random 
potential can be thought of as being smooth over a dis- 
tance of the width L. With periodic boundary conditions 
in x-direction one would again conclude that eigenstates 
are plane waves in x-direction, labeled by N different 
center coordinates Y q (q is an integer times 2ir/L) in y- 
direction, and the eigenvalues E q would be determined 
by the value of the random potential at center coordi- 
nate Y q . This scenario is also consistent with Eq. JT^), 
because the Fourier transform of the random potential at 
V(Y q ) yields the matrix elements h m (see also Eq. ||). In 
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V. CONCLUSIONS 
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FIG. 11. Squared amplitudes of a typical state in the ran- 
dom band matrix model for situation (B) with strongest cor- 
relation between matrix elements. The state is represented in 
the two-dimensional Landau representation. 

the absence of periodic boundary conditions the situation 
changes. For energies far from the band center one ex- 
pects, that the corresponding eigenstates are localized on 
equipotential contours of the random potential and are 
centered at some value Y q . However, close to the energy 
band center eigenstates become extended and one typical 
eigenstate is shown in Fig. 11. Although this state has 
a preferred orientation in x-direction it is by no means 
localized to a small region in y-direction. It fluctuates 
strongly, it has non- vanishing values all over the system, 
and it also shows large areas of low probability. There- 
fore, the multifractal exponent ao is larger than in the 
standard quantum Hall situation. 

Let us try to give heuristic arguments of how to estimate 
the value of ao- For that purpose we recall that, quite 
generally £qid is of the order of the number of transverse 
modes N c times the relevant scattering length I (for a 
discussion see e.g. [Q). In our situation N c = N and 
I w d ~ L. Therefore, the quasi-lD localization length 

5/3 

is estimated to be ~ L y , and it is much larger than L y 
p3[ . We may thus assume that the state is critical and 
has, in the strip-representation, a value ao ~ 2.26 when 
the fractal analysis is restricted to sizes much larger than 
d ~ L. Recall that we have chosen the Landau represen- 
tation corresponding to an aspect ratio a = 1. Therefore, 
the value of ao found by box counting in that represen- 
tation must be different. The box counting method uses 
squares of size I 2 in the 2D Landau representation with 
0=1. This corresponds to rectangular boxes in the strip- 
representation, where the length in y direction scales as 
the third power of the length in x-direction. Thus, the 
"effective volume" is Z 4 . By this reasoning ao{C = N) — 2 
will be, in the a = 1 representation, two times larger than 
in the strip-representation, and we may expect that we 
should find a C ~ N ^ « 2.54 by box-counting in 2D Lan- 
dau representation with a = 1. Indeed, this estimate is 
compatible with our findings as displayed in Fig. 10. 



We have studied a novel type of matrix models, the cor- 
related random band matrices. We used numerical di- 
agonalization and performed a multifractal analysis to 
analyze the localization-delocalization properties of such 
matrix models in the thermodynamic limit of infinite ma- 
trix size. 

The parameters of correlated band matrices are the band 
width B and the correlation parameter C. We offered 
three interpretations: (i) independent quantum parti- 
cles on a one dimensional chain with correlated hop- 
ping, (ii) independent quantum particles on a quasi-one- 
dimensional strip with correlated coupling of channels, 
and (iii) in some range of its parameters the models 
resemble two dimensional quantum Hall systems. For 
B ~ C ~ y/W a transition from localized to critical states 
in the band center occurs, and the corresponding critical 
exponents are close to those of real quantum Hall sys- 
tems. Furthermore, we found the following qualitative 
behavior when keeping the band width ~ V7V constant: 
A reduction of correlations suppresses multifractality (i.e. 
criticality) at the band center and finally, for C = 1, the 
ordinary non-critical random band matrix ensemble is 
reached which shows localization lengths £ ~ B 2 . In- 
creasing correlations beyond C ~ VN, the transition to 
critical states in the band center remains, however their 
multifractality seems to be more pronounced. The frac- 
tal critical exponent for extreme correlations, C(N) = N, 
turned out to be compatible with a heuristic estimate. 
Therefore, our numerical results suggest that the corre- 
lated band matrix models show transitions from local- 
ization to critical derealization on approaching the en- 
ergy band center, provided the band width scales like 
B(N) ~ y/~N and the correlation parameter scales like 
C(N) ~ N* with 1/2 < t < 1. It should be pointed 
out that correlations lead to stronger localization off the 
band center, while they lead to critical derealization at 
the band center. 

We hope that our work initiates more studies on the en- 
semble of correlated random band matrices with a general 
behavior of B(N) ~ N s , C(N) ~ where s, t may vary 
between and 1, and to reach solid statements about the 
localization behavior in the thermodynamic limit. We 
also like to point out, that the "standard quantum Hall 
case" of the correlated random band matrix models is 
not only a simple matrix realization for quantum Hall 
systems, but has a very interesting distinction from other 
representative models for the quantum Hall universality 
class (for an overview over such models see |3^]). The 
correlated random band matrix does not incorporate any 
handedness related to the magnetic field. This handed- 
ness is essential in all other representative models that al- 
low for the existence of extended states. In the correlated 
random band matrix model, however, the connection to 
a quantum Hall system goes via the Landau representa- 
tion, which takes the handedness into account. Fortu- 
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nately, the question of localization and derealization is 
not restricted to that representation. In the correlated 
band matrix model the correlation of matrix elements is 
the key for the quantum Hall transition. It would be very 
interesting to construct a manageable field theoretic for- 
mulation for the correlated random band matrix model. 
This may be possible when taking advantage of the fact 
that the correlations are given by constraints which may 
be included by Lagrangian multipliers. 
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